Magnetic excitations in underdoped Ba(Fei_ x Co x ) 2 As2 with x =0.047 
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The magnetic excitations in the paramagnetic-tetragonal phase of underdoped 
Ba(Feo.953Coo.o47)2As2, as measured by inelastic neutron scattering, can be well described 
by a phenomenological model with purely diffusive spin dynamics. At low energies, the spec- 
trum around the magnetic ordering vector Qafm consists of a single peak with elliptical shape 
in momentum space. At high energies, this inelastic peak is split into two peaks across the 
direction perpendicular to Qafm- We use our fittings to argue that such a splitting is not due 
to incommensurability or propagating spin- wave excitations, but is rather a consequence of the 
anisotropics in the Landau damping and in the magnetic correlation length, both of which are 
allowed by the tetragonal symmetry of the system. We also measure the magnetic spectrum deep 
inside the magnetically-ordered phase, and find that it is remarkably similar to the spectrum of the 
paramagnetic phase, revealing the strongly overdamped character of the magnetic excitations. 



I. INTRODUCTION 

The discovery of high-temperature superconductiv- 
ity in iron-based compounds^ triggered an intense re- 
search activity to understand these materials 3 4 . By now, 
there is ample evidence indicating the unconventional 
nature of the superconducting state, such as the exis- 
tence of a magnetic resonance mode^ the observation 
of half-integer flux-quantum in niobium-iron pnictide 
loops^ the microscopic coexistence between magnetism 
and superconductivity and the linear temperature- 
dependence of the low-T superfluid density in some 
compounds^. Due to the proximity of the superconduct- 
ing transition temperature (T c ) to a magnetic instabil- 
ity, one of the main candidates to explain this unconven- 
tional state is the pairing mechanism mediated by spin 
fluctuations^. 

In order to investigate the viability of spin fluctuations 
as the pairing mechanism, it is paramount to have a con- 
sistent and complete description of the magnetic excita- 
tions. Inelastic neutron scattering (INS) measurements 
performed by several groups have provided a common 
picture for the so-called 122 compounds (^4Fe2As2, A de- 
noting alkaline earth)M^UH Similar to the cuprates, the 
iron pnictide superconductors are characterized by high 
magnetic energy scales extending beyond 200 meV. Such 
high energy scales could possibly support high tempera- 
ture superconductivity where electrons are paired by spin 
fluctuations.^ In the iron pnictides, the striped mag- 
netic ground state is characterized by the ordering vector 
Qafm = (|, |,1) (in tetragonal notation), with antifer- 
romagnetic spin correlations in the direction parallel to 
(|, |, 0) and ferromagnetic spin correlations along the di- 
rection perpendicular to (§, §,0). 



Unlike the Neel (or checkerboard) AFM ordered state, 
this striped AFM state breaks the tetragonal symme- 
try of the lattice, and results in an orthorhombic dis- 
tortion with the longer (antiferromagnetic) axis a paral- 
lel to Qafm and the shorter (ferromagnetic) axis b per- 
pendicular to Qafm- However, in several systems, the 
tetragonal symmetry-breaking, probed by the develop- 
ment of the orthorhombic distortion, takes place above 
the magnetic transition temperature In, implying that 
long-range magnetic order is not the origin of the struc- 
tural transition. Rather, it is possible that magnetic fluc- 
tuations can spontaneously break the tetragonal symme- 
try in the paramagnetic (PM) phase , giving rise to the so- 
called electronic nematic stateP^L^ Thus, the systematic 
study of the magnetic spectrum via INS is important not 
only to shed light on the superconducting pairing mech- 
anism, but also to clarify the nature of the tetragonal 
symmetry-breaking in the iron pnictides. 

INS studies have focused on the magnetic excita- 
tion spectrum in both the parent compounds, which 
are antiferromagnetic (AFM) metals, and doped su- 
perconducting (SC) compositions that are character- 
ized by strong paramagnetic spin fluctuations. In each 
case, the INS measurements of the spin excitations 
have revealed surprisingly strong in-plane anisotropics. 
For instance, the orthorhombic-AFM state of the par- 
ent CaFe2As2 compound displays a significant in-plane 
anisotropy in the spin- wave velocities and in the zone- 
boundary energies^^. When the spin wave spectrum 
is fitted with the J\ — J<2 Heisenberg model, the only 
way to account for such anisotropics is to use very dif- 
ferent values for the nearest-neighbor exchange interac- 
tions in the directions parallel (Ji a ) and perpendicular 
(Jib) to Qafm - to the extent that one is negative (anti- 
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ferromagnetic) and the other positive (ferromagnetic) 
While Ji a ^ Ji6 is justified in the orthorhombic-AFM 
phase, the disparity between their values, especially con- 
sidering the rather small orthorhombic distortion a — 6, 
suggests that a simple description in terms of localized 
moments may not suffice. Even in the tetragonal-PM 
phase of CaFe2As2, where the low-energy spin fluctua- 
tions were shown to be diffusive in nature, the magnetic 
correlation lengths perpendicular and parallel to Qafm 
are unequal. Such an in-plane anisotropy does not con- 
tradict the tetragonal symmetry of the high-temperature 
phase as Qafm = (§, §,1) is a point of two- fold sym- 
metry located on the Brillouin zone boundary. Similar 
anisotropic spin fluctuations have also been observed in 
the tetragonal phase of the Ba(Fei_ ;E Co :E )2As2 samples 
with x = 0.065 15 and 0.074 12 , where long -range AFM or- 
der is absent. Whether the observed anisotropics of the 
spin-wave velocities in the orthorhombic-AFM phase and 
of the magnetic correlation lengths in the tetragonal-PM 
phase share a common origin remains an open issue. 

At high energies (> 80meV), anisotropic spin fluctu- 
ations within the Fe layer develop a transverse splitting 
in paramagnetic Ba(Fei_ :E Co :E )2As2 with x = 0.074.^ 
This unusual spectral feature has also been recently ob- 
served in the parent BaFe2As2 compound Although 
different groups generally agree on these observations, 
the nature of the high-energy peak-splitting is currently 
under intense debate. The work by Li et a/., which 
first observed such splitting, attributed it to the exis- 
tence of quasi-propagating modes with different velocities 
along the transverse and longitudinal directions^!, with- 
out, however, providing an explanation for their highly 
anisotropic character. Based on band-structure calcu- 
lations, Park et al. associate the high-energy splitting 
to an incommensuration of the SDW ordering vector^, 
whose detection apparently falls below the resolution 
limit in the low-energy regime. Harriger et al. per- 
formed a detailed study of the high-energy spectrum in 
both the orthorhombic and tetragonal phases of the par- 
ent compound BaFe2As2p^. Starting from an effective 
overdamped anisotropic J\ a — Jit, — J2 model, they ar- 
gue that the high-energy cross section cannot be fitted 
with Ji a = Jib and isotropic Landau damping even in 
the tetragonal PM phase. As a result, they suggest that 
the tetragonal symmetry is broken at high energies and 
that the peak-splitting is a signature of nematic short 
range order. A similar detailed analysis was made by 
Ewings et al. on another parent compound SrFe2As2pM. 
They find that while an effective Ji a — J15 — J 2 model 
with fixed parameters is unable to describe the data for 
different temperatures and energies, a band-structure cal- 
culation similar to Park et a I™ correctly describes the 
complete spectrum of magnetic excitations, without re- 
quiring tetragonal symmetry-breaking. 

In this paper, we study the spin fluctuation spectrum 
at high energies in Ba(Fei_ ;E Co :E )2As2 with x = 0.047 us- 
ing INS. In combination with the results in the literature 
on x = 0, 0.065, and 0.074, the new data allows for a 



discussion of the systematic evolution of the high energy 
spin susceptibility as a function of both composition and 
temperature covering undoped, underdoped, and opti- 
mally doped samples. We find that the transverse split- 
ting is relatively insensitive to composition, suggesting 
that it is impervious to fine tuning of the Fermi surface 
and unlikely to arise from static incommensurability^. 
The splitting is also insensitive to temperature, including 
temperatures above and below Tn, i.e., it is insensitive to 
the presence of long-range AFM order. This temperature 
independence was used to suggest that the spin fluctua- 
tions are characteristic of the orthorhombic phase at all 
temperatures, implying a tetragonal symmetry-breaking 
above T/v 14 . Here, we adopt the opposite point of view, 
and suggest that the transverse splitting is a consequence 
of the anisotropic spin fluctuations already present in 
the tetragonal-PM phase, manifested by anisotropics in 
both the magnetic correlation length and Landau damp- 
ing. A similar conclusion was drawn by Refs. [15] and [21] 
using first-principle calculations combined with random 
phase approximation (RPA) and dynamical mean field 
theory (DMFT) approaches, respectively. Here, we sup- 
port our claims by developing a phenomenological low- 
energy model for the dynamical magnetic susceptibility 
derived from an effective three-band model for the Fe 
sublattice. The main advantage of this approach is to 
provide a simple description of the anisotropic spin fluc- 
tuations in terms of a few physically relevant parameters. 
We are able to fit all our INS data using a single set of 
parameters, showing the suitability of this phenomeno- 
logical model as a description of magnetic spectrum of 
the iron pnictides. 

Our paper is organized as follows: details of our exper- 
iment are contained in Sectio n \H\ observed spin excita- 
tions are described in Section |III| effects of composition 
are discussed in Section \LV\ the models used to fit the 
observed spin excitations are detailed in Section |V| de- 
tails of our analysis of the neutron scattering data are 
presented in Section [Vl} finally, our results are discussed 
in Section [ylll 



II. EXPERIMENTAL 

Inelastic neutron scattering measurements were per- 
formed on the ARCS 22 spectrometer at the Spallation 
Neutron Source at Oak Ridge National Laboratory on 
Ba(Fei_ ;E Co :E )2As2 with x = 0.047. The sample con- 
sists of ten co-aligned single-crystals grown from excess 
FeAs and Co As, as outlined in Ref. 23, with a total 
mass of ~ 2 grams mounted with a horizontal [HHL] 
scattering plane (in tetragonal notation). The sam- 
ple undergoes a tetragonal to orthorhombic structural 
transition at T$ = 60 K, orders antiferromagnetically 
at Tn = 47 K and magnetic long-range order coexists 
with superconductivity below T c = 17 K. The sample 
was aligned with the c-axis along the incident beam di- 
rection and measurements were performed with an in- 
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cident neutron energy of Ei = 250 meV. The inelastic 
scattering spectra were measured at T = 5K (i.e. in 
the orthorhombic, AFM and superconducting state) and 
T = 70 K (i.e. in the paramagnetic, tetragonal state). 
Throughout the manuscript, the neutron scattering data 
are described in the tetragonal IA/mmm coordinate sys- 
tem with Q = ^ (H + K) i + ^ (H - K) j + ^Lk = 
(H + K,H-K, L) where a = 3.95 A and c = 12.95 A. In 
tetragonal I4/mmm notation, Qafm = (| 5 [^ = ^? 
if=0]. if and K are defined to conveniently describe 
diagonal cuts in the IA/mmm basal plane as varying H 
(K) corresponds to a longitudinal [H, H] scan (transverse 
[K,—K] scan) through Qafm- 



III. DESCRIPTION OF THE SPIN 
EXCITATIONS IN THE x = 0.047 SAMPLE 

The elastic magnetic scattering and low-energy spin 
dynamics (< 15meV) of this same sample were mea- 
sured by triple-axis neut ron scattering techniques and 
published elsewhere! 24 * 25 1 The x = 0.047 sample has a 
small ordered moment (~ 0.2/ie) and the low-energy spin 
dynamics in the AFM phase at T = 25 K (normal state, 
T c < T < Tn) are characterized by strong damping. Un- 
like similar measurements on the ^4Fe2As2 parent com- 
pounds where a sharp spin gap appears in the ordered 
state, the spin gap is strongly broadened or absent for 
the x = 0.047 composition, suggesting that the damping 
energy scale is comparable to the spin gap energy in the 
ordered state (~ lOmeV). Time-of-flight measurements 
are required in order to ascertain if well-defined collec- 
tive spin wave excitations appear in this system at higher 
energy transfers. 

Time-of-flight measurements of the spin excitations at 
T = 5 K are shown in Fig. [I] The magnetic excitations 
are observed to emanate from Qafm and are steeply dis- 
persive, extending to energies approaching 150 meV. Sur- 
prisingly, we were not able to observe sharp and well de- 
fined collective spin wave excitations at T = 5K at any 
energy despite the long range AFM order. This observa- 
tion confirms that strong damping is a general feature of 
the spectr um, in agreement with the low-energy triple- 
axis resultsP^EU The dampi ng of t he spin waves is larger 
than that found in CaFe2 As J 19 | 2Q | and BaFe2As2^ where 
long-lived collective modes are seen above ~ 50 meV and 
provide clear evidence of a conical spin wave dispersion. 
In addition a well-defined spin gap is also observed in 
CaFe2As2 and BaFe2As2- In this respect, the spin fluctu- 
ation spectrum for x = 0.047 is more like that measured 
in the optimally doped and paramagnetic compositions 
of Ba(Fei_ a; Co :c )2As2 (x > 0.06), where the spin fluc- 
tuations are short-ranged and diffusive in character ! 12 * 26 * 

The diffusive nature of the spin fluctuations for x = 
0.047 in the AFM state is further emphasized by mea- 
surements in the paramagnetic state at T = 70 K > In 
(Ts), shown in Fig. [5] which display a magnetic spectrum 
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FIG. 1: (color online) Inelastic neutron scattering data mea- 
sured on Ba(Feo.953Coo.o47)2As2 at T — 5K using ARCS 
with E{ — 250 meV and the crystal aligned with the inci- 
dent neutron beam along the [0,0,1] direction, (a) Trans- 
verse slice of the data along the [K, —K] direction through 
Qafm = (§, \,L) after averaging over H = 0.5 ± 0.05 r.l.u. 
(b) Longitudinal slice of the data along the [if, H] direction 
through Qafm = (|, after averaging over K — =b 0.1 

r.l.u. In each panel, the color scale represents the intensity of 
scattered neutrons. 



nearly identical to the one at T = 5 K. Thus, we find 
that despite long-range AFM order, the response func- 
tion has a diffusive character, with the spin fluctuations 
being characterized by strong damping and short correla- 
tion length as opposed to sharp collective spin-wave exci- 
tations. In addition, a comparison of the transverse and 
longitudinal cuts shown in Figs. [2|a / )-(d / ) and (a // )-(d // ) 
confirm a substantial anisotropy of the in-plane correla- 
tion lengths and a transverse splitting of the high energy 
spin fluctuations both above and below Tn- Similar fea- 
tures were reported previously for the x = 0.074 com- 
position. The similarity of the spin fluctuations at the 
two temperatures allows us to average the data taken 
at T = 5K and 70 K in order to improve the counting 
statistics for subsequent analysis described below. Fig- 
ure [3ja) and (b) show temperature-averaged spectrum 
cut in the transverse and longitudinal directions, respec- 
tively, while panels (c)-(h) show constant energy slices of 
the data in the (if, K) plane. In particular, the slices at 
80 (Fig.^e)) and 100 meV (Fig.|Ji)) clearly display the 
transverse splitting of the spin excitations. 
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H= 0.5zb0.06 r.l.u. 



K= 0±0.15 r.l.u. 
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FIG. 2: Comparison of inelastic neutron scattering data mea- 
sured on Ba(Feo.953Coo.o47)2As2 at T = 5 K (empty circles) 
and T — 70 K (filled circles). Transverse cuts of the data at 
(a) 100 ± 10 meV, (b') 80 ± 10 meV, (c') 60 ± 10 meV, (d') 
40 ± 10 meV. Longitudinal cuts of the data at (a") 100 ± 
10 meV, (b") 80 ± 10 meV, (c") 60 ± 10 meV, (d") 40 ± 10 
meV. 



IV. COMPOSITION DEPENDENCE 

The high-energy spin fluctuations have been stud- 
ied in several compositions of the Ba(Fei_ x Co^ se- 
ries (x = (P, 0.06^1 0.07^, 0.0^51, and the x = 
0.047 composition discuss ed here) as well as in the par- 
ent compound CaFe2 As2 ! n l 19 l 2Q I In each composition, 
the spin excitations are seen to extend up to energies 
> 100 — 150 meV. In the paramagnetic phase, the spin 
fluctuations in momentum space display a pronounced 
elliptical anisotropy within the Fe layer, with softer mag- 
netic excitations extending along the [K, —K] -direction 
(transverse to Qafm) and stiffer excitations along the 
[H, ^-direction (longitudinal to Qafm)- In this regime, 
the anisotropy is defined in terms of the magnetic corre- 
lation lengths along antiferromagnetic (longitudinal cor- 
relation length, £l) an d ferromagnetic bonds (transverse 
correlation length, £t) of the stripe magnetic structure: 
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(i) 



- corresponds to isotropic spin fluctuations while 
rj£ = 1 (rj£ = —1) corresponds to the extreme limit 
where £l ^> £t (£l <^ £t)- This two-fold anisotropy 



at Qafm = is allowed by the simple tetrago- 

nal symmetry of the Fe sublattice. 13 The origin of the 
anisotropy can be deduced from both localized^ and itin- 
erant descriptions of the magnetism in these materials. 
In the limit of local-moment magnetism, one can use the 
J\ — J2 model to show that rj£ — J\ / 2 J2^ Within this 
approach, 77^ = corresponds to no coupling between 
the two interpenetrating Neel sublattices making up the 
stripe AFM structure (Ji = 0), whereas r/^ = 1 corre- 
sponds to the classical stability limit of the stripe ordered 
antiferromagnetic state, which gives way to G-type mag- 
netic order for J\ > 2J2. In the itinerant approach, rj£ is 
a consequence of the ellipticity of the electron pockets. 

For those compositions having long-range AFM order, 
the spin-waves dispersion is anisotropic within the Fe 
layer. To compare the anisotropy in the AFM and PM 
states, we define the anisotropy of the spin wave velocity 



as 



Vc 



(2) 



where cl (ct) is the spin wave velocity along the longi- 
tudinal (transverse) direction, respectively. If one uses 
a local- moment description for the spin wave dispersion, 
z.e., the Ji a — Jib — J2 model, and if the orthorhombic 
distortion is small (ao ~ fro), then r] c is given by 

r] — — — (3) 
4J2 + Jia — Jib 

When Ji a « J^ — Ji, the spin- wave velocity 
anisotropy reduces to the same result as that obtained 
for the correlation lengths in the tetragonal paramagnetic 
phase, T] c = 77£ = Ji/2J 2 . 

Using these expressions, we can compute experimental 
values for the correlation length and spin wave veloc- 
ity anisotropics for several different iron arsenide com- 
positions in both the ordered and paramagnetic phases, 
as shown in Table [TJ In those systems with long-range 
AFM order, the anisotropy of the low-energy spin fluctu- 
ations does not change strongly above Tn, i.e., between 
the AFM-orthorhombic and PM-tetragonal phases. This 
observation casts much doubt on the "nematic spin fluid" 
model proposed in Ref. [14j In that model, local or- 
thorhombic distortions are proposed to exist based on the 
analysis of short-wavelength zone boundary spin-waves 
at Q = (1,0, L) using the Ji a — Jib — J2 model. The 
temperature independence of these zone boundary spin- 
waves suggests that J\ a and J15, which are very different 
in the orthorhombic phase, remain so even in the tetrag- 
onal paramagnetic phase. However, this local symmetry 
breaking cannot hold in the long- wavelength limit, where 
the crystallographic symmetry must be respected and 
an average nearest-neighbor exchange 2Ji = Ji a + Jib 
must be assumed. Such a scenario would predict a siz- 
able change in the momentum space anisotropy of low 
energy spin excitations above and below Tn, since 

" tCt = 1 + Jla Z Jlb (4) 
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FIG. 3: (color online) Inelastic neutron scattering data measured on Ba(Feo.953Coo.o47)2As2 using ARCS with E- x — 250 meV 
and the crystal aligned with the incident neutron beam along the [0,0,1] direction. Data at T = 5K and 70 K are summed 
together to improve statistics, (a) Transverse slice of the data along the [K, —K] direction through Qafm = (|, after 
averaging over H — 0.5 ±0.05 r.l.u. (b) Longitudinal slice of the data along the [H, H] direction through Qafm = (§, §, L) after 
averaging over K — =b 0.1 r.l.u. Constant energy slices in the (if, fT)-plane averaged over an energy range of (c) 35 =b 5meV, 
(d) 50 ± lOmeV, (e) 80 ± lOmeV, (f) 100 ± lOmeV, (g) 120 ± lOmeV, and (h) 150 ± lOmeV. In each panel, the color scale 
represents the intensity of scattered neutrons, where the maximum intensity in each panel, in arb. units, is (a) 10, (b) 10, (c) 
10, (d) 8, (e) 5, (f) 3, (g) 2, (h) 1; and the minimum intensity in all panels is -1 arb. units. 



Using the values of the exchange constants J\ a and 
Jib extracted from the spin- wave fitting of BaFe2As2, we 
have ^ et /^° rt ~ 1.7. Thus, we would expect a change 
of the in-plane anisotropy by 70%, however, no change 
of the anisotropy is observed experimentally. 14 This tem- 
perature independence suggests that the J\ — J2 model is 
an appropriate parameterization of the spin dynamics in 
the long- wavelength limit (in either the tetragonal or or- 
thorhombic phases). At shorter wavelengths/higher ener- 
gies, the local-moment models clearly run into trouble^ 
and a description that explicitly takes into account the 
itinerancy of the magnetic moments is needed.^ It is 
interesting to note that spin fluctuations become more 
anisotropic with Co doping in the Bal22 system, reach- 
ing values 77^ > \ for x > 0.065. 

At energies > 80 — 100 meV, the anisotropic spin fluc- 
tuations split along the transverse direction and seem 
to form separate counter-propagating branches. We note 
that this splitting is unusual and not typical of propagat- 
ing spin waves, where constant energy contours will form 



a ring of scattering from collective modes that propagate 
in all directions. This transverse splitting phenomenon 
was first noted in the x = 0.074 composition^ and has 
since been observed in x = parent compound. 14 In the 
current study on x = 0.047, we also observe the trans- 
verse splitting, as shown in Figures [2] and [3j We note that 
the energy where this transverse splitting takes place and 
its intensity are temperature-independent (Fig. i.e., 
the splitting is unmodified even deep in the AFM ordered 
state, a feature also observed in BaFe2As2™ Viewed as a 
propagating mode where ct = huj/ Aq, Table[l]shows that 
the velocity of this mode is weakly dependent on compo- 
sition within the Bal22 series, with a value of 316 meV A 
at x = that softens to ~ 240meVA for optimally 
doped superconducting compositions x > 0.065. The 
weak dependence of this splitting on composition casts 
some doubt on its interpretation as an incipient incom- 
mensurability, i.e., a band nesting effect, which should 
be very sensitive to composition!^ Recently, static in- 
commensurate AFM order was indeed observed, but 
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TABLE I: Composition dependence of the transverse prop- 
agation velocity and anisotropy of spin fluctuations in 
A(Fei_ cc Co cc )2As2 (A = Ca, Sr, Ba). For undoped compounds 
in the antiferromagnetic ordered (AF) phase, the transverse 
velocity and anisotropy parameters are obtained from a fit 
to a model of propagating spin waves. For paramagnetic (P) 
compounds, the parameters are obtained from either damped 
spin wave or diffusive models. See references for further de- 
tails. 
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200 


PM 
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only in a very narrow range of Co concentration from 
0.056 < x < 0.06P 



V. MODELING THE SPIN EXCITATIONS IN 
THE PARAMAGNETIC PHASE 

A. Diffusive model 

To describe the magnetic spectrum in the paramag- 
netic state, we first consider an itinerant model that cap- 
tures the diffusive character of the spin dynamics. One 
can perform a first-principles calculation for the spin dy- 
namics by using the complete density functional theory 
(DFT) band struct ure a nd incorporating the electronic 
interaction via 

RPA H3ii5J or DMFTP Indeed, such calcu- 
lations provide a good description of the magnetic spec- 
trum of different compounds, but their complexity makes 
it difficult to sort out the essential physics responsible for 
the behavior of the magnetic spectrum. 

Here, instead of considering the full band structure in- 
put, we adopt a phenomenological approach. The low- 
energy magnetic spectrum of this model has been dis- 
cussed in Refs. 7 and 29. Near Qafm m the lA/mmm 
tetragonal coordinate system, we obtain the complex dy- 
namical susceptibility: 



x(q + Qafm,^) = 



Xo 



(5) 



Here, £ is the bare magnetic correlation length and xo 
is the inverse of the typical magnetic energy scale. The 
imaginary term appears due to the decay of spin excita- 
tions (paramagnons) in particle- hole pairs, and 7 q is the 



corresponding momentum-dependent Landau damping. 
The momentum-space anisotropy term rj^ is an indirect 
consequence of the ellipticity of the electron pockets con- 
sidered in this model. Taking the imaginary part, which 
is the quantity measured by INS, we can cast the result 
in the form: 



x" (q + Qafm,^) = 



xo£ 2 ^r q 



■r*(i 



with 



and 



F[l + a 2 (q 2 + 2rirq x q y )] 



(6) 



(7) 



(8) 



In Eqn. (8), we expand the momentum-dependent 
damping T q for small deviations from Qafm, introduc- 
ing the anisotropy parameter 7?r- To ensure that T" 1 is 
maximum at q = 0, we restrict our analysis to \rjr\ < 1. 
Similarly, to ensure that the magnetic correlation length 
in Eqn. ^ is well-defined, we assume |%l< 1- Notice 
that Eqn. ^ naturally gives rise to Eqn. ([T]) considered 
before. 

At low energies, this form of the susceptibility cor- 
rectly captures the elliptical shape of \" (q> ^) peaked 
at Q = Qafm(<1 = 0), in agreement with the INS 
measurements on a variety of 122 compounds! 11 * 12 * The 
sign of the ellipticity depends on 77^; for the parent and 
Co-doped BaFe2As2 compounds, such as the one shown 
in Fig. [3j the longitudinal correlation length [parallel 
to Qafm = (§,§)] £l = £\A + % is longer than the 
transverse correlation length [perpendicular to Qafm = 
(k->\)\ £t = £\A — implying that rj£ > 0. On the 
other hand, the sign of 77^ depends on the band structure 
(see, for instance Ref. Il3|) , and there is a priori no reason 
for it to be positive or negative. This is to be contrasted 
to the classical AFM J1 — J2 model, which always predicts 
77 oc J1/J2 > 0. 

At high energies, as uo increases, the momentum- 
dependence of the Landau damping eventually leads to 
the splitting of the single peak of \" (q?^)- Along the 
transverse direction, the low-energy peak splits into two 
symmetric peaks for energies uj > cjt, with: 




2£ 2 (l - gg) 
a 2 (1 - Vr) 



(9) 



Similarly, along the longitudinal direction, the peak- 
splitting takes place for energies above 




2g 2 (l + ^) 
a 2 (l+ Vr )- 



(10) 
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K(r.l.u.) H(r.l.u.) H (r.l.u.) H(r.l.u.) 

FIG. 4: (color online) Calculations of the inelastic neutron scattering spectrum as obtained from a model of itinerant diffusive 
spin fluctuations, as described by Eqn. §6§. The parameters of the model were obtained by fits to the neutron scattering data 
for Ba(Feo.953Coo.o47)2As2 as shown in Figs. [6] and [7] (a) Transverse slice of the data along the [K,—K] direction through 
Qafm = (§, §, L) after averaging over H = 0.5 ±0.05 r.l.u. (b) Longitudinal slice of the data along the [H, H] direction through 
Qafm = (§, \ ,L) after averaging over K — ± 0.1 r.l.u. Constant energy slices in the (H, K)-p\ane averaged over an energy 
range of (c) 35 ± 5 meV, (d) 50 ± 10 meV, (e) 80 ± 10 meV, (f) 100 ± 10 meV, (g) 120 ± 10 meV, and (h) 150 ± 10 meV. In each 
panel, the color scale represents the intensity of scattered neutrons, where the maximum intensity in each panel, in arb. units, 
is (a) 10, (b) 10, (c) 10, (d) 8, (e) 5, (f) 3, (g) 2, (h) 1; and the minimum intensity in all panels is 0. 



Therefore, if rj£ = rjr, the elliptical shape of x" (q?^) 
at low energies evolves in an elliptical-ring structure for 
uj > co>t = co>l- On the other hand, for 7^ > rjr (rj£ < 
rjr) the splitting happens at a lower energy along the 
transverse (longitudinal) direction. Thus, depending on 
the magnitude of the difference Au = uo^ — ^l, there can 
be a wide regime of energies where the low-energy single- 
peak response of x" (q? ^0 a t Qafm splits into two peaks 
equidistant from Qafm- For Auo > (i.e. rj£ > rjr), these 
two peaks split along the direction transverse to Qafm, 
whereas for Auj < (i.e. rj^ < rjr), the two peaks split 
along the longitudinal direction to Qafm- Therefore, the 
experimental data in the parent and Co-doped BaFe2As2 
samples (Fig. |3| are compatible with r]£ > rjr- 

Figure [4] shows calculations of the neutron scattering 
cross-section for the diffusive model using parameters de- 
termined from fits of the neutron data for the x = 0.047 
compound to Eqn. The fits to the data are described 
in detail below. This figure shows the same slices of the 
neutron intensity as in Fig. [3] and can be compared di- 



rectly to the data. In general, we find that the neutron 
data is best described with anisotropy parameters rj£ > 
and rjr < 0, resulting in anisotropic ellipsoids of scatter- 
ing at low energies and a transverse splitting at higher 
energies. 

B. Ballistic model with propagating modes 

Besides the purely diffusive dynamics discussed in the 
previous section, one can consider a different phenomeno- 
logical form for the spin susceptibility. In the so-called 
ballistic model, the dynamics is governed by propagating 
overdamped spin-waves in the disordered paramagnetic 
state 

x"(q + Q A FM,w) = x ° eujT £2 (li) 

q 

with anisotropic spin- wave velocity: 
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FIG. 5: (color online) Calculations of t he inelastic neutron scattering spectrum as obtained from a ballistic model of overdamped 
spin waves, as described by Eqn. (11). The parameters of the model were obtained by fits to the neutron scattering data for 
Ba(Feo.953Coo.o47)2As2 as shown in Figs. [6] and [7] (a) Transverse slice of the data along the [K, —K] direction through 
Qafm = (§, \ ,L) after averaging over H = 0.5 =b 0.05 r.l.u. (b) Longitudinal slice of the data along the [H,H] direction 
through Qafm = (§, §, L) after averaging over K — =b 0.1 r.l.u. Constant energy slices in the (H, i^)-plane averaged over an 
energy range of (c) 35±5meV, (d) 50±10meV, (e) 80±10meV, (f) 100±10meV, (g) 120=bl0meV, and (h) 150±10meV. In 
each panel, the color scale represents the intensity of scattered neutrons, where the maximum intensity in each panel, in arb. 
units, is (a) 10, (b) 10, (c) 10, (d) 8, (e) 5, (f) 3, (g) 2, (h) 1; and the minimum intensity in all panels is 0. 



whereas here we assume a linear dispersion. 



2 2 2/2 

c^q = c (q 



^r] c q x q y ) 



(12) 



The anisotropy in the velocity, r] Cl is the same as that 
introduced in Eqn. (J2|. Presumably, the spin wave veloc- 
ities in the paramagnetic phase are comparable to those 
found deep in the AFM ordered state. The form of the 
susceptibility in Eqn. (11) is diffusive in nature (consist- 
ing of relaxational dynamics and a single peak response) 
at low energies, where the spin- wave wavelength is longer 
than the correlation length {i.e. uj < c q £ _1 ). The renor- 
malized spin-wave modes will appear in the form of a 
broad elliptical ring of scattering at constant uo when 
^ > c q£ -1 is satisfied. The high-energy form of the 
scattering then represents a section of the damped, con- 
ical spin-wave dispersion. This approach was used to 
describe the anisotropic quasi-propagating mode postu- 
lated in Ref. 12. A similar approach was also used in 
Ref. [14j although in that case a local moment Heisenberg 
model was employed to describe the spin wave dispersion 



Figure [5] shows calculations of the neutron scattering 
cross-section for the ballistic model using parameters de- 
termined from fits of the neutron data for the x = 0.047 
compound to Eqn. (11). The fits to the data are de- 



scribed in detail below. This figure shows the same slices 
of the neutron intensity as the data in Fig. [3J as well as 
the diffusive model in Fig. |4j We see that the anisotropy 
of the spin wave velocity with r] c > describes the el- 
liptical shape of the scattering cross-section. At high 
energies, the scattering assumes the expected form of an 
elliptical ring. We note that while the spectra in pan- 
els (a) and (b) of Figs. [2jj5] appear similar, there is a 
qualitative difference in the higher energy slices in pan- 
els (c)-(h). The diffusive model correctly reproduces the 
transverse splitting observed by experiment, whereas the 
ballistic model can only produce elliptical rings or, in the 
extreme case of r] c « 1, an elliptical ring that is pinched 
in the middle (see Fig. |5jg) and (h)). 
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VI. ANALYSIS OF THE NEUTRON 
SCATTERING DATA 

We now describe the procedure used to fit the INS 
data to both the diffusive and ballistic models. The in- 
elastic neutron scattering data is directly related to the 
imaginary part of the magnetic susceptibility through the 
relation 



5(Q,w) oc / 2 (Q)x"(Q,w)(l - e-fc-ATj-i (13) 

Each scattered neutron collected at ARCS is recorded 
as an individual event, where the time- of- flight and de- 
tector position is stored. It is impossible to visualize the 
data without first converting its format from a list of 
flight times and pixel positions to a spectrum for each 
detector pixel, this necessitates collecting the data into 
energy bins such that each spectrum is a histogram of 
the number of detected neutrons per unit energy trans- 
fer. This process is often referred to as data reduction 
and was in this case accomplished with the DANSE soft- 
ware packaged 

The observed inelastic neutron scattering intensity 
7(Q,co>) includes contributions from various non-sample 
sources - primarily the crystal holder and the sample 
environment. We also consider non-magnetic scattering 
from the sample itself as a source of background. The 
collective intensity from these sources comprise the back- 
ground function B(Q,lj). The background function was 
estimated by averaging intensity from equal Q and u sec- 
tions of the detector which were far from magnetic inten- 
sity. 



Y,j B(Q j ,uj j )5(Qi - Qj)S(uJi - ujj) 



(14) 



where the sum goes over all detectors far from mag- 
netic intensity and the Qij^ij correspond to bin cen- 
ters. This approximate to B(Q,uj) can only properly 
account for contributions to the background which are 
isotropic about the incident beam direction and is, more 
accurately, B(Q,uj). An estimate of the sample-scattered 
intensity can be obtained by subtracting the averaged 
background estimate from the observed intensity. 



S(Q,u) = I(Q,u) - B(Q,u) 



(15) 



The background subtracted spectra for detectors 
within the range 0.25 < H < 0.75r.l.u., -0.5 < K < 
0.45r.l.u. and 20 < E < 150 meV were collectively fit 
with the Levenberg-Marquardt algorithm to each model 
function, Eqns. (|6|) and JTT1). In addition to either of the 



model functions, a background was fit to account for an 
observed residual component of £?(Q,co>) caused by a de- 
tector bank top-to-bottom asymmetry of the measured 
intensity, the fit residual background function was 



where the ai are positive model parameters and K (in 
r.l.u.) is a measure of the vertical displacement from the 
horizontal detector plane. Our fitting routine ensured 
that, for all cases, |^| < 1. Volumetric H,K,E data was 
fit in order to avoid artifacts introduced by binning which 
are inherent in 1-D cuts and 2-D slices. Fit parameter 
values and their associated errors, plus values of 77^ de- 
rived from the fit parameters are presented in Table [TTJ 
We observed that fits of the data were rather insensitive 
to changes in r]r , and subsequently fixed r]r to — 1 for our 
fitting. This limit gives ^-independent (constant) damp- 
ing in the longitudinal direction (r) and ^-dependent 
damping in the transverse direction (T q = T(l-\-2a 2 q 2 )). 



Diffusive (Eqn. 


6) 


Ballistic (Eqn. 11 


) 


Xo£ 2 (arb. units) 
T (meV) 

a (A) 

£t (A) 

a (A) 
rjr 


5.0(4) 
10.7(11) 
9.4(5) 
5.0(3) 
2.0(2) 
-1 


Xo£ 2 (arb. units) 
T (meV) 

a (A) 

£t (A) 

c(meVA) 

Vc 


6.3(7) 
8.5(11) 
11.8(9) 

8.2(6) 
450(40) 
0.73(6) 




0.57(6) 




0.35(9) 



Br(Q,w) 



-CLq 



(16) 



TABLE II: Fit parameter values and calculate d 7fe for the two 
models: diffusive (Eqn. |6j), and ballistic (Eqn. [Tlj). 



This manuscript describes inelastic neutron scattering 
data from Ba(Feo.953Coo.o47)2As2 which exhibit a trans- 
verse splitting in the spin excitations which has been 
previously observed for other doping concentrations of 
Co. A detailed understanding of the spin fluctuations 
is an important component of understanding the ori- 
gin of superconductivity in the doped BaFe2As2 com- 
pounds. We find that this high energy splitting is present 
in the presence of both a paramagnetic-tetragonal and 
antiferromagnetic-orthorhombic-superconducting phase. 
Furthermore, we have successfully described our data by 
means of a model involving anisotropic Landau damping 
and magnetic correlation length, which does not break 
the high temperature tetragonal symmetry of the system. 
Our findings show that it is not necessary to introduce 
symmetry breaking, quasi-propagating modes, or incom- 
mensurability in order to explain the high energy split- 
ting of the spin excitations in Ba(Fei_ ;E Co :E )2As2, which 
readers of PRB should find important and interesting. 

We use Eqn. ( p~3| ) for the neutron intensity to test both 
the diffusive and the ballistic models, Eqns. (|6| and ( pTj) 
respectively. The fittings are presented in Figs. [6] and [7| 
(red curve for the diffusive and blue curve for the bal- 
listic model) for fixed energy and fixed momentum, re- 
spectively. For energies below 60meV, the two models 
give nearly identical line shapes, which is not surprising, 
since both models reduce to an identical description of 
diffusive spin excitations at low energies. On the other 
hand, at higher energies, the purely diffusive model is in 
better agreement with the data, since the longitudinal 
cut has a maximum at q = whereas the transverse cut 
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H= 0.5±0.06 r.l.u. 



K= 0±0. 15 r.l.u. 




<><> 



K= 0.00±0. 025 r.l.u. 
H= 0.5±0.3 r.l.u. 



FIG. 6: (color online) Fitting results for the diffusive (red 
lines) and ballistic (blue lines) models to the summed T = 
5K and 70 K data (open symbols). Transverse cuts (') and 
longitudinal cuts (") through Q = (§, §, L) at (a) E = 120 ± 
lOmeV, (b) E = 100 ± lOmeV, (c) E = 80 ± lOmeV, (d) 
E 55±5meV, (e) E 45±5meV, and (f) E 35±5meV. 



has a minimum at q = 0. The ballistic model, on the 
other hand, gives the opposite: a minimum at q = for 
the longitudinal cut and a maximum at q = for the 
transverse cut. 



VII. DISCUSSION 

According to the fits of the INS data, neither the purely 
diffusive nor the ballistic model is strongly favored over 
the other. Mostly, this is due to the broad and over- 
damped nature of the spin excitations and the conver- 
gence of the two models at low energies. Surprisingly, 
even though both models are low-energy descriptions of 
the magnetic excitation spectrum, the fitting to the data 
reveals that they are able to capture some of the intricate 
physics responsible for the high-energy behavior. Physi- 



o- 

5-b ' 




1 1 I 1 1 1 1 I 



K= 0.05±0.025 r.l.u. 
H=0.5±0.3 r.l.u. " 




100 150 
E (meV) 

FIG. 7: (color online) Fitting results for the diffusive (red 
lines) and ballistic (blue lines) models to the summed T = 
5 K and 70 K data (open symbols) . Spectra cuts centered 
at (a) (0.5, 0.5, L), (b) (0.55, 0.45, L), (c) (0.6, 0.4, L), (d) 
(0.65, 0.35, L), and (e) (0.7, 0.3, L). 



cally, the two models represent two different conceptual 
pictures for AFM fluctuations in the paramagnetic phase. 
The ballistic model would represent spin wave propaga- 
tion in disordered AFM, whereas the diffusive model in- 
cluding Landau damping is based more closely on a quasi- 
particle description obtainable from either the simplified 
3-band structure discussed here, or by first principles cal- 
culations in Refs. 13 and T3] However, the qualitative 
form of the scattering throughout the (#, K) -plane, in 
particular the observation of the split transverse modes 
at high energies, seems to favor a description of the exci- 
tations in terms of the diffusive model, without the need 
to introduce propagating modes (compare Figs. |2j [3] and 
[4|. In the diffusive model, the peak splitting is caused by 
the interplay between the momentum-space anisotropics 
associated with the Landau damping and the magnetic 
correlation length, both of which follow naturally from 
the band structure of the iron pnictides. 

It is interesting to note that a similar discussion regard- 
ing the existence or not of propagating modes took place 
also in the context of a much simpler material, namely, 
bcc iron (see, for instance, Refs. [3D and 32). In that 
case, claims for the existence of propagating modes at 
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low energies followed from observations of splitting in 
the constant-energy cuts of the INS data similar to our 
case. It was subsequently shown that the nature of the 
damping in an itinerant ferromagnet, where T q — >• as 
q —> 0, can result in a splitting that resembles a spin- 
wave dispersion.^ In the ferromagnetic case, energy cuts 
at constant-momentum should show an inelastic peak in 
the ballistic model, whereas no such peak will be observed 
in the diffusive model. In our specific case of paramag- 
netism in the stripe AFM ordered phase, one can see 
(Fig. [7]) that both the ballistic and purely diffusive mod- 
els produce an inelastic peak in the constant-momentum 
energy scan. Thus, unlike the ferromagnetic case, the 
presence or absence of an inelastic peak in the spectrum 
cannot discern between the two models. 

The excitations in the paramagnetic phase of the par- 
ent BaFe2As2 compound have been described using a 
J\ — J2 Heisenberg model combined with a phenomeno- 
logical momentum-dependent damping function. This 
approach was able to describe the spin fluctuations close 
to the zone boundary only after allowing for tetrago- 
nal symmetry-breaking, as discussed in Ref. [T4J As de- 
scribed above, this "nematic spin fluid" model is called 
into question since the proposed model should lead to 
large changes in the low energy anisotropy of the spin 
fluctuations below Tn, which are not observed. In or- 
der to describe the peculiar transverse splitting, the au- 
thors chose a strongly anisotropic form for the Landau 
damping that resulted in large damping along the lon- 
gitudinal direction and small damping in the transverse 
direction. This choice effectively washes out the longitu- 
dinal spin waves, leading to the split transverse modes. 
In our notation, the authors chose rjr < — 1 meaning that 
the damping function has a saddle point around Qafm, 
rather than a minimum, as one would expect. While sim- 
ilar in spirit to the approach considered here, our model 
does not require local symmetry breaking in the param- 
agnetic state. We note, however, that our analysis does 
not preclude the existence of a nematic phase in the iron 
pnictides. To probe the nematic phase, it would be more 
appropriate to perform INS measurements in detwinned 
samples (see, for instance [33] and also l34|) . 

Finally, we comment on the role of incommensuration 
of the SDW ordering vector. In Ref. [T3l first princi- 
ples calculations suggest that an incommensurability is 
present for the parent compound BaFe2As2- It was then 
proposed that the high-energy peak splitting would be 
nothing but a manifestation of this incommensurability. 
If indeed the magnetism in the iron pnictides is of itin- 



erant nature, then it is reasonable to expect, on general 
grounds, the development of incommensurability at some 
critical doping (see the seminal work of Rice^). Subse- 
quent to those predictions, static incommensurability was 
indeed observed in a narrow range of doping concentra- 
tions in the Co-doped Bal22 system. 28 The incommensu- 
rability is observed to develop in a first-order fashion with 
doping concentration) 2 ^ However, in the INS data pre- 
sented here for the underdoped Ba(Feo. 953C00. 047)2 As2, 
as well as in the previously presented data for the 
parentPS and optimally doped Ba(Feo.926 Coo.o74)2As2p', 
the high-energy splitting of the inelastic peaks occurs at 
roughly the same energy u sp n t « 80 meV and momentum 
in both the commensurate AFM ordered state and the 
PM state. Thus, it unlikely that the transverse splitting 
is associated with any instability towards incommensu- 
rate order. 

In summary, we developed a simple expression able 
to capture the main properties of the spectrum of the 
magnetic excitations in these materials. In particular, 
we showed that the observed transverse splitting of the 
inelastic peaks of \" (q> u ) a ^ high energies, as well as 
their dispersions, can be attributed to the interplay be- 
tween the anisotropics of the magnetic correlation length 
and of the Landau damping. This explanation follows 
naturally from the model band structure we used, and 
does not require breaking of the tetragonal symmetry, 
or the introduction of quasi-propagating modes or in- 
commensurability. We demonstrated the applicability 
of this model to the 122 iron pnictides by fitting low 
and high energy INS data on the underdoped compound 
Ba(Fe .953Co . 047)2 As 2 . 
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